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The equations for the electromagnetic two-body problem are neutral-delay equations that for 
generic initial data have solutions with discontinuous derivatives. If one wants to use these neutral- 
delay equations with arbitrary initial data, solutions with discontinuous derivatives must be allowed. 
Surprisingly, this same neutrality is compatible with the recently developed variational method 
with mixed-type boundaries for the Wheeler-Feynman electrodynamics. We show that two-body 
electromagnetic orbits with discontinuous velocities are physically necessary by showing that orbits 
with vanishing far-fields amost everywhere must have some discontinuous velocities on a few points. 
We generalize the Wheeler-Feynman electrodynamics with the variational method to include all 
continuous trajectories, allowing piecewise-differentiable weak solutions represented by trajectories 
with fields defined almost everywhere (but on a set of points of zero measure where velocities jump). 
Along with this generalization we formulate the generalized absorber hypothesis that the far-fields 
vanish asymptotically almost everywhere and show that bounded two-body orbits satisfying the 
generalized absorber hypothesis need to have discontinuous derivatives on a few points. We also give 
the general solution for the family of bounded non-radiating two-body orbits. We discuss the physics 
of orbits with discontinuous derivatives and show that these conserve the physical momentum, 
stressing the differences to classical variational methods. Last, we discuss how the electromagnetic 
variational method with mixed-type boundaries is well-posed but lacks reversibilty. 
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INTRODUCTION 



Since the speed of light is constant in inertial 
frames, the equations of motion for point-charges are 
state-dependent-delay equations. More specifically, the 
Wheeler-Feynman electrodynamics [l| has mixed- type- 
state-dependent-neutral-delay equations of motion for 
the two-body problem. The theory of delay equations 
is still on the make 0-01 but it is already known that 
delay-only equations with generic continuous initial data 
have continuous solutions with a discontinuous derivative 
at the initial time. The derivative becomes continuous 
at the next breaking point Q and progresses from C k to 
C fc+1 at successive breaking points. On the other hand, 
for a neutral-delay equation with generic continuous ini- 
tial data @, 0] solutions can have discontinuous deriva- 
tives at all breaking points. If one wants to use the elec- 
tromagnetic neutral-delay equations with arbitrary ini- 
tial data, solutions with discontinuous derivatives must 
be expected. Surprisingly, this same neutrality is com- 
patible with the recently developed variational method 
with mixed-type boundaries for the Wheeler-Feynman 
electrodynamics Two-body electromagnetic orbits 
with discontinuous velocities are not that odd physically, 
and in fact here we show that any sequence of orbits radi- 
ating less and less converges to orbits with discontinuous 
velocities. For these limiting orbits where the accelera- 
tion is not defined on a few points, the variational method 
*5j] offers a well-posed alternative to define weak trajecto- 



ries beyond those satisfying a neutral-delay equation ev- 
erywhere. The variational method || describes naturally 
all continuous orbits with possibly discontinuous deriva- 
tives. The relation to Maxwell's electrodynamics is of in- 
clusion, i.e., this generalization contains the C 2 orbits 
of the Wheeler-Feynman theory. As shown in Ref.[f|, if 
boundary data are such that the extremum orbit turns 
out to be piecewise C 2 with continuous velocities, the 
equations of the Wheeler-Feynman electrodynamics hold 
everywhere with the exception of the breaking points (a 
set of measure zero). Here we show that if the velocity 
is discontinuous at these breaking points, further non- 
local momentum conservations are required by the vari- 
ational method, another striking difference to the clas- 
sical local two-body variational methods. Finally, if the 
extremum in not piecewise C 2 the variational method 
yields a generalized weak dynamics beyond description 
by the Wheeler-Feynman neutral-delay equations. 

The variational equations of Ref. [H along piece- 
wise C 2 orbits include the electromagnetic-fields as the 
farthest-reaching couplings to the other particle's trajec- 
tory. In the original articles, Wheeler and Feynman [l[ 
attempted to derive an electrodynamics with retarded- 
only fields from the hypothesis that the universal far- 
fields vanish at all times [H, yielding a self- force depend- 
ing on the derivative of the acceleration. Here we gen- 
eralize the absorber hypothesis [l[ to trajectories with a 
discontinuous velocity and acceleration on a set of mea- 
sure zero by defining the far-fields almost everywhere 
with the usual electromagnetic formulas while allowing 
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the fields to be undefined on a set of points of measure 
zero, thus arriving at the generalized absorber hypothesis 
that the far-fields vanish almost everywhere, henceforth 
called G.A.H. The generalized flux of the Poynting vec- 
tor is a physical property involving an integral, so that 
it can be evaluated even with fields undefined on a set of 
measure zero. Here we show that a bounded two-body 
G.A.H. orbit must have discontinuous derivatives on a 
few points, by giving the general solution to a simple 
neutral-delay equation. 

We henceforth adopt a unit system where the speed of 
light is c = 1 and apply the usual formulas of electro- 
dynamics to piecewise-defined trajectories with the ex- 
ception of the few points where the past/future veloc- 
ities/accelerations are undefined, where we declare the 
fields are undefined (a set of measure zero for the flux 
integral). The far-electric field of a point charge in the 
Wheeler- Fey nman electrodynamics involves a sum of ad- 
vanced and retarded fields, 



while the far-magnetic field is given by 



B = ~n+ x E a 

2 + 



dv 



-n_ x E ret , 
2 



(1) 



(2) 



where unit vectors n± point away from the ad- 
vanced/retarded position of the charge, respectively [fj. 
We consider a spatially bounded orbit enclosed by a 
sphere of radius R in an inertial frame, with R much 
larger than the maximum orbital excursion, so that 
n + = n_ = n. The space-time points (i, Rn) on the 
sphere are specified by the time t and the unit vector n 
and the Poynting vector P — E x B evaluated with Eqs. 
(HJ) and © at (t, Rn) becomes 
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adv 1 2 



IE 



ret 1 2 



}n 



(3) 



where single bars denote Euclidean modulus, and we have 
used the transversality of the far-fields, i.e., n-E re * = 
n . E adi, _ q Notice that the G.A.H. implies the vanish- 
ing of the flux © because |E ret | 2 = |E a(to | 2 = almost 
everywhere, while the vanishing of the surface integral of 
Eq. <j3j> alone does not imply the G.A.H. For example the 
circular two-body orbits [3 H| of the Wheeler- Fey nman 
electrodynamics do not satisfy the G.A.H. We henceforth 
introduce an index k = 1, 2 to label the charges. The far- 
magnetic field is proportional to ~E r k et (t, n) by 



Br(*,n)=nxE£ e 



(4) 



For a continuous and piecewise C 1 trajectory x,t : 
[— Afc, T^f] — > K 3 the retarded far-electric field of charge 
qk defined piecewise at the space-time point {t, Rn) is 
given by the Lienard-Wiechert formula [6J 



K et (t,n)=q k 



n x {[n - v fc _] x a fc _} 
(1 - n ■ v k .) 3 R ' 



(5) 



Because of Eq. ((3]) it suffices to study the vanishing of 
the electric far-retarded fields and we henceforth assume 
the orbit is time-reversible so that the vanishing of the 
retarded far-fields implies the vanishing of the advanced 
far- fields. In Eq. ([5]), unit vector n points from the 
charge's retarded position x/^ifc-)) to the space- 
time point (t, Rn) while v,t_ = dx.k/dtk\t k _ and = 
d 2 Xk/dtl\ tk _ are respectively the Cartesian velocity and 
Cartesian acceleration of the point charge, which along 
a generalized continuous-only orbits are defined almost 
everywhere in past light-cone of (t,Rn). The time of 
particle k in lightcone with (t, Rn) is defined implicitly 
by the retardation condition 



t k -=t - |x fe (£fe_) - Rn\, 



(G) 



where single bars stand for Cartesian distance. If the 
piecewise defined velocity Vfc_ = dx.k/dtk\ t - is lesser 
than light, Eq. (fTT|) defines the retarded time tk- as 
an implicit function of time t with a piecewise defined 
derivative 
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(7) 



dt (l-n-Vfe_) 
We can use Eq. J7|) to express the far-electric-field (|5|) as 



Er(*,n) = ^x^[x fc (t fc _)], 



(8) 



where Xfc(ifc_) is the position of particle k at time t k -. 
Henceforth charge 1 is supposed positive and equal to q 
while charge 2 is negative and equal to — q. The G.A.H. 
for piecewise diffcrcntiable orbits is expressed almost ev- 
erywhere by 



E 



E^ 



Er=^x^(x 1 (t 1 _)-x 2 (t 2 _))=0. 

(9) 

Assuming the Xi(ii) piecewise C 2 , condition ((9} has a 
piecewise-linear continuous solution defined for t in each 
closed-open interval t e (to-,^o-+i] with creZ by 

XiOi_) - x 2 (t 2 _) = D CT (n) + nf a (t, n) + (t - t a )V a {n), 

(10) 

where the D CT (n) and V^n) are arbitrary bounded func- 
tions and the f a (t, n) are bounded and piecewise C 2 . It 
is possible to choose n • D cr (n)= and adjust D CT (n) to 
make the left-hand-side of Eq. (jTU)) continuous. 

Along a spatially bounded orbit Eq. (j6|) is approxi- 
mated for large values of R by 



t k -=t - R + n ■ x fe (t fc _) 



(11) 
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Notice that Eqs. (jlip yield an implicit relation between 
ti and 4J, 

ii_ - i 2 - = n- (x 1 (t 1 _) - x 2 (4 2 _)). (12) 

It is instructive to use the retardation conditions ([6]) with 
k = 1 , 2 and the implicit function theorem to express 4 2 as 
a function of ti and n. We define the influence interval of 
point (i 2 — , x 2 (f 2 _)) by the interval containing 4i_ when 
n varies arbitrarily in Equation (|12[) . i.e., 

4 2 _ - |xi(ii) - x 2 (4 2 _)| < h < t 2 - + |xi(ti) - x a (f a -)|. 

(13) 

The time span (fl3|) is from the retarded light-cone 
time of (4 2 _, x 2 (4 2 _)) to the advanced light-cone time 
of (4 2 _,x 2 (4 2 _)Y, as with the mixed- type boundary con- 
ditions of Ref . [5J . Notice that the future light-cone 
appeared naturally in the two-particle problem, even 
though we were dealing only with the retardation con- 
ditions [flT]). It follows from Eqs. and {jTUJ) that 

/ CT (4,n) = (ti_ -t 2 _) - (t-On-Vitn). (14) 

and we can re-write Eq. (flT))) as 

xi(*i_)-x 2 (* 2 -) = D„(n)+((i_-t 2 _)n-((-t ff )nx(nxy c 

(15) 

Since linear growth is unbounded, the only globally 
C 2 orbit must have V a (n) = Vcr and it follows from Eq. 
([T5]l with V CT (n) = that 



xi(ti_) - x 2 (< 2 _) = D CT (n) + - i a _)n. 

The derivative of Eq. (fll))) respect to time yields 
vi_ v 2 _ 



where 



(l-n-vi_) (l-n-v 2 _) 



(l-n-vi_) (l-n-v 2 _) 



(16) 



(17) 



(18) 



Notice that Kyi is still given by Eq. (jT8|) and with 
nonzero V a (n) the right-hand-side of Eq. (|20|) forms a 
complete 3-dimensional basis to express any vector (in- 
side or outside the plane of Vi(4i_) and v 2 (4 2 _) ). Equa- 
tion ([12")) still allows one to move n in a cone with axis 
along Xi(t]_) — x 2 (4 2 _) 7^ in a way that fixes 4i_ and 
4 2 _ while t changes with Eqs. (fTTj) . By choosing ifi 2 
and a nonzero V^-(n) for each t G (4^, f^+i] we can de- 
scribe any vector on the left-hand-side of Eq. (|20l) . so 
that there is no inconsistency. Example time-reversible 
orbits satisfying Eq. (121)1) are piecewise-constant-velocity 
orbits generated by jumping one velocity at a given time 
while the other velocity jumps either in the backward or 
forward light-cone times symmetrically, as well as at ev- 
ery time in the forward and backward light-cones of a 
discontinuity time (the sewing chain illustrated in Fig. 
2 of Ref. [5j ). These piecewise-linear polygonal orbits 
can be checked to satisfy Eq. © by direct substitution 
and use of Eq. (|7J) . To find the most general solution of 
Eq. (|15l) we notice that for any given piecewise defined 
trajectory x 2 (4 2 _) and given D CT (n) and V^n) Eq. (fT5|) 
determines in general only a function xi(ti_,n) of the 
two variables (4i_,n) by 
( n ))- 

xi(ii_,n) = x 2 (t 2 _ )+D ff (n)+(ti- -« 2 - )n- )nx (nx (n)). 

(21) 

The implicit function theorem further determines 4 2 _ and 
t as functions of 4i_ and n by Eqs. (fTT|) . (fl"2|) and l[2T]) . 
A physical trajectory must satisfy the extra consistency 
requirement that the Xi(4]_,n) determined by Eq. (|l?l"j) is 
a function of 4i_ only, i.e., 



9xi(fi_,n) 
dn 



0. 



(22) 



Equation (|12|) allow us to move n in a cone with axis 
along xi(ti_) — x 2 (t 2 _) 7^ in a way that fixes t\- and 
4 2 _ while changing t with Eqs. (fTTj) . On the other hand, 
for fixed 4x_ an d *2- the left- hand-side of Eq. (flTf spans 
a plane of the fixed vectors Vi_ and v 2 _, so that Eq. 
(TT7)) ) can hold only if K 12 — 0, which combined with 
Eqs. ([17} and JTHJ) yields 



Condition (|22l) applied to the right-hand side of Eq. (p?Tj) 
is the extra condition determining a consistent trajec- 
tory. Since (|2"2")) must hold for all values of 47/ inside each 
interval of the piecewise defined orbit, we must also have 
inside each piecewise interval that 

d 2 xi(ti„,n 



dti-dn 



d dxi(4i_,n) 
<9n l flti_ j 



0. 



(23) 



which can be expressed as 

d 9t 2 _(4i_,n) 

— (v 2 _-n) 

an a4i_ 



<94(4i_,n) 
-n ^- -nx(nx^(n)] = 0. 

Oil— 



vi(4i_) = v 2 (4 2 _). 



(19) 



Equation (fT9"| defines piecewise-constant velocities, and 
the only bounded choice is Vi = v 2 = 0, as discussed in 
[loj . Nontrivial alternatives necessitate the introduction 
of a few discontinuities, i.e., use Eq. (H"5j) with V a ^ 0, in 
which case the piecewise derivative of Eq. (fT5|) respect 
to time yields 



(24) 

A symmetric compatibility condition follows by exchang- 
ing 1 and 2 in Eq. (|24[) . determining a set of equation 
whose general solution is 

54 2 _(4i_,n) <94(4i_,n) 
[v 2 _ - nj £ + n — -n x (nxK,(n) 



dti- 



dti- 



Vl_ 



v 2 - 



(l-n-vi_) (l-n-v 2 _ 



JsTi 2 n - n x (nxV/ 7 (n)). 

(20) 



A CT (4i_), 

54i_(4 2 _,n) <94(4 2 _,n) 
[vi_ - n] ^ + n H ^ -n x (nxV^(n) 



dt 2 



dU 



B a {h-), 



(25) 
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for arbitrary functions and B G (t 2 _). We stress 

that the discontinuities introduced by the V a (n) play 
an essential part in the solution; Trying to solve ei- 
ther one of Eqs. ([231) with V a (n) — yields the for- 
mer glitch generated by the rotation of n with fixed t\- 
and i2- outside either the plane of A a {tiJ) and V2- or 
the plane of B a (t 2 _) and vi_. Otherwise with nonzero 
V a (n) there is a complete basis to endow rotational in- 
variance to Eqs. (|25|) . which define the most general non- 
radiating bounded orbit. It can be seen that piecewise- 
linear polygonal orbit are a special case for constant A a 
and B a . Next in line would be to find the extremum 
bounded orbits of the variational method [f| that also 
satisfy Eqs.(|25|). which should give an even better ap- 
proximation for the orbits studied in Ref.[t| (even more 
astonishing hit at the spectroscopic lines!). The phys- 
ical need of trajectories with discontinuous velocities is 
the justified as limitind orbits defiuned by sequences of 
bounded two-body orbits radiating less and less. 

Now that we know that bounded G.A.H. orbits must 
involve a discontinuous derivative on a few points, as 
compatible also with neutral-delay equations, let us ex- 
amine the variational method of [5|. We henceforth as- 
sume that Let a continuous trajectoris Xi(ti) and x 2 (t 2 ) 
S R 3 and piecewise C 1 in the above defined (t a ,t a+ i] 
intervals. The mixed-type boundary conditions for the 
variatiuonal method [B| as illustrated in Figure 1 are the 
initial point Oa for trajectory 1 plus the segment of tra- 
jectory 2 inside the lightcone of Oa, and the endpoint 
Lg for the trajectory 2 plus the segment of trajectory of 
particle 1 inside the lightcone of Lb- For variations of 
trajectory 1 the variational method has a LagrangianQ 



Si = 



£i(xi,vi,x 2 , v 2 )dii 



(26) 



-miy 1 — vfdti 

(1 - vi ■ v 2+ ) 
2ri2+(l - ni2+-v 2+ ) 

(1 - Vl • V 2 _) 

2ri2_(l - ni2--v 2 _) 



where unit vector ni2± points respectively from the ad- 
vanced/retarded position (t 2 ±, x^fei)) ^° ^ e s P ace_ time 
point (ii,xi(ii)) and v 2 ± = dx 2 /dt 2 \ t2± while r 12+ = 
xi(ii) — x 2 (^2±)| are the distances in light-cone to tra- 
jectory 2. The variational method is defined in the 
space of C 1 orbital variations of trajectory 1 with fixed 
endpoints, i.e., 



<5xi(ii = 0) =0, 
5xi (h = T L -)=0. 



(27) 



The linear expansion of Si about the orbit defines the 
Frechet derivative, i.e., 



SS i = / Kir 1 ■ <* x i) + • <* v i)]^i + o(|^i| 2 ), 

(28) 

where (a -b) is the scalar product of R 3 . In particular 
when x\ : [0, Tj,-] — > R 3 is piecewise C 2 we can integrate 
Eq. ()28[) by parts in each interval and re-arrange yielding 

SSx = r\^-[^-i^Wt. (29) 



1 , chci dt 9vi 

(T—N — 1 „ r 



The extremal orbits in the class of continuous orbital 
variations must (i) Satisfy the Euler-Lagrange equation 
to vanish the first line of Eq. (|29[) with arbitrary varia- 
tions and (ii) Conserve the momentum currents 



dLi.t. 
<9 Vl h ° 



dLi 
dvi 



(xi, Vl 



dLi 
dvi 



x x ,v la ) = 0. (30) 



For local Lagrangians usually one has dL/dvi = 
G(xi,v 1 .x 2 ) such that Eq. (f3T))) requires the velocity 
to be continuous. For piecewise C 2 orbits continuity of 
velocity implies that only globally C 2 are possible, a pop- 
ular result of classical mechanics. On the contrary, for 
the Lagrangian given in Eq. (14) of Ref [B| the current is 



dLi 
9vi 



miVi 



(31) 



V2- 



V24 



2ri2-(l - ni2--v 2 -) 2ri2+(l - n 12+ -v 2 ^ 



with an analogous continuity along the trajectory of par- 
ticle 2. The Noether two-body-momentum of Ref. 
formula (A23), involves an integral that is insensitive to 
velocity jumps plus two non-local currents given by Eq. 
(|3T|). (see Eqs. (A25) and (A26) of Ref. (3) that are sen- 
sitive to jumps. Therefore even the physical Noether two- 
body momentum is conserved as long as Eq. pip is con- 
tinuous across the jumps. The physical justification to 
generalize electrodynamics with the variational method 
of Ref. J5[ is precisely to include the non-radiating orbits 
of Ref.[9j in the family of bounded G.A.H. orbits. We 
have here sinthesized three different reasons for orbits 
with a discontinues derivative : (i) inclusion of bounded 
G.A.H. orbits, (ii) the conservation of Noether 's momen- 
tum (iii) compatibility with the neutrality of the equa- 
tions and (iv) the fact that the variational method is nat- 
ural in a space completed to contain orbits with discon- 
tinuous velocities. The linearization of the neutral-delay 
equations of motion used in Ref-Q about circular orbits 
can produce piecewise-defined orbits with discontinuous 
derivatives. Moreover, these orbits can be approximated 
by a series of stiff approximants having continuous deriva- 
tives and an increasingly sharp edges to approximate the 
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velocity jumps, i.e., the stiff perturbations in Ref.Q that 
predicted orbits in the atomic magnitude with a surpris- 
ing precision from a non-radiation condition involving a 
flux integral. We notice that this surprising difference 
to the popular variational principles of classical mechan- 
ics required a minimum of two bodies and a non-local 
Lagrangian. 

Electromagnetism was originally formulated with the 
integral laws of Ampere, Gauss and Faraday, and only 
much later differential equations holding everywhere were 
introduced by Maxwell. The requirement of a second 
derivative existing everywhere is not needed for particle 
dynamics, where one is concerned only with an integral of 
the force along trajectories. The variational method with 
mixed-type boundaries of the Wheeler-Feynman electro- 
dynamics [f| allows for such generalization, a step back 
from Maxwell's equations in the sense of weak solutions. 

The bounded two-body G.A.H. orbit is a seed to start 
a perturbation theory for a three-body motion by ex- 
tending the variational method of [f| to three charges 
and placing the third charge at a large distance. For 
the three-body variational method the interaction terms 
fall sufficiently fast (i.e.. the far- fields of the bounded or- 
bit) almost everywhere along the third trajectory. We 
see that a bounded G.A.H. orbit is special in the sense 
that it does not "disturb" the universe. Moreover, there 
is an stability provided by the G.A.H. when the flux of 
the semi-sum field vanishes in a large universe contain- 
ing many charges. The G.A.H. is much better than hav- 
ing the (non-zero) fluxes of both retarded and advanced 
fields cancel each other. If the G.A.H. holds, an offending 
perturbation of size e in the retarded far-field (or in the 
advanced far-field) perturbs the flux Eq. ^ at 0(e 2 ), 
otherwise the flux ([3]) is perturbed at 0(e). 

The dependence on the mixed-type boundaries must 
be investigated for scattering trajectories with discon- 
tinuous velocities and accelerations at the boundaries; 
These are likely to have future continuations involving 
stiffcr jumps at later times, so that particles collide with 
laboratory boundaries, which can be regarded as a gen- 
eralized type of radiative loss. This is to be contrasted 
with the fate of C°° trajectories determined by continu- 
ation of C°° initial data as discussed in Ref. [Hj|:- For 



these it is shown in [11] that the equations of motion with 
advanced fields can be solved for the most-advanced ac- 
celeration of the other particle, yielding delay-only equa- 
tions, Eq. (22) of Ref. [12j. so that no future information 
is needed to define solutions. On the contrary, the future 
data in general is only part of the mixed-type boundary 
data used in [B| . The variational method Q with generic 
mixed-type boundary data determines orbits with discon- 



tinuous velocities in a way that backward continuation is 
not unique. The lack of well-posed reversibility solves 
a long-standing debate for the absorber hypothesis with 
C°° trajectories; The generalization of Ref. [IH allows 
no unique backward continuation because there are no 
equations of motion to run backwards. To continue the 
trajectories into the past with the variational method 
one needs to guess pre-historical boundaries for one par- 
ticle and solve another mixed-type-boundary- value prob- 
lem with that postulated segment of trajectory plus the 
other particle's earlier past data. This construction is 
not unique because there is not a unique pre-history. The 
lack of a unique past for the other particle does not vio- 
late causality if the predicted future of one particle is far 
enough from the past data of the other particle because 
the advanced/retarded fields involve respectively a future 
or a past position measured at another spatial point by 
another clock synchronized a la Einstein. In the theory 
of relativity one is not allowed to compare times not mea- 
sured at the same point so that to falsify the predicted 
future of the other spatial point one would have to travel 
to that spatial point. For mixed- type boundary data Q 
with a large enough time-separation between past data 
and future data, any piecewise-subluminal orbit arrives 
after the predicted future happens, so that no contra- 
diction is involved, only ill-posed (non-unique) backward 
continuation. 
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